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Wind profiles in and over trees 
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Abstract: One of the most important and frequently studied variable in forests and the most basic element in governing transport 
processes of airflow is wind speed. The study of wind profile, defined as the change of wind velocity with height, and wind ve¬ 
locity are important because of tree physiological and developmental responses. Generally, wind profiles above the ground or at 
a canopy surface follow classical logarithm law, but wind profiles in a single tree and in a forest stand are not logarithmic. This 
paper summarizes the results of wind profile studies within a single tree, in a forest stand, above the forest canopy and in a 
forest area from recent research in a coastal pine forest. The results demonstrate that: 1) wind profiles with in a single conifer 
tree crown showed an exponential function with height, 2) wind profiles in forest stands were able to be expressed by attenuation 
coefficient of wind, 3) wind profiles over a forest canopy could be determined using profile parameters (friction velocity, rough¬ 
ness length and displacement), and 4) for a forest area, the extreme wind speed could be predicted reasonably using the 
methods developed for the design of buildings. More research will be required to demonstrate: 1) relationships between wind 
profiles and tree or stand characteristics, 2) the simple methods for predicting wind profile parameters, and 3) the applications of 
wind profile in studies of tree physiology, forest ecology and management, and the detail ecological effects of wind on tree 
growth. 
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Introduction 

One of the most important and most frequently studied 
variables in forests and those most basic in governing 
transport processes and dispersion of airflow are wind 
speed and temperature (Raynor 1971; Zhu et al. 2000). 
Wind causes extensive forest damage and large amounts 
of soil being wind eroded and transported to another loca¬ 
tion, and it is also a major factor influencing other micro¬ 
climate factors such as the heat and vapor exchange near 
the soil surface (Coutts and Grace 1995; Peltola et al. 
2000), the transfer processes of insects and disease, and 
the growth of trees (Hannah et al. 1995). In all of these 
interactions, the wind profile, i.e., the change of wind ve¬ 
locity with height, is very important, because both physio¬ 
logical and developmental responses of trees to wind are 
sensitive to the shape of the wind magnitude and height. 
Therefore, much research has been carried out to under¬ 
stand the physical processes involved so as to enhance 
silvicultural practices and forest management (Landsberg 
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and James 1971; Oliver 1971; Kotoda and Hayashi 1980; 
Cionco 1985; Ishizaki and Ota 1987; Gardiner 1995; Liu & 
Kotoda 1996; Lemone etal. 1998; Zhu etal. 2004). 

The variation of wind speeds with height from a surface 
can be expressed empirically by various formulae. Gener¬ 
ally, for ground level studies and for moderate or strong 
winds, the wind profile expressed by equation (1) is valid 
(Raynor 1971; Oliver and Mayhead 1974; Lyles and Allison 
1979; Kimmins 1987; Abtew et al. 1989; Dyrbye and Han¬ 
sen 1997; Maki 1999; Zhu etal. 2001a). 

U z =u,~ ln(^^.) + <5, h) 

K Z 0 

where 14 is the horizontal mean wind velocity (m-s' 1 ) at 
height z, u* is the friction velocity (m-s' 1 ), K is von Kar- 
man’s constant (0.40, dimensionless); z is the height of 
measurement from a reference plane (m); do is the 
zero-plane displacement (m); z 0 is the roughness parame¬ 
ter (m); and <D Z indicates that mean velocity also depends 
on Richardson's number (Ri, dimensionless). 

The wind profile (equation 1) has been applied in many 
studies. When we consider the motion of air in a tree or a 
forest stand, it becomes quite apparent that airflow within 
the tree crown and forest canopy is a significant aspect of 
the surface boundary layer (Cionco 1985). 

The aerodynamic characteristics of forest stand have 
been an interest to forestry and agriculture for a long time, 
particularly in regard to the classic problems such as wind 
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damage, seed transport, disease and insect control, and 
fire prevention. More complex approaches to explore the 
atmospheric constituents that are affected by airflow inter¬ 
acting with vegetative structures have also been carried out. 
Several meteorological models have been constructed and 
much data collected (Reifsnyder 1955; Taketa 1965; Allen 
1968; Raynor 1971; Kotoda and Hayashi 1980; Amiro 
1990a, b; Lemone et at. 1998). It has thus been shown that 
the wind profile is the most sensitive indicator of the degree 
and effectiveness of turbulence (Bull and Reynolds 1968; 
Gardiner 1995). Additionally, wind profile is also important 
in assessing wind damage or stand stability (Blackburn and 
Petty 1988b; Peltola 1996b; Gardiner et al. 1997), and in 
examining wind effects on microclimates in forest stand, 
e.g., on the distribution of rainfall, light and some gases 
(Telewski 1995; Kuraji et al. 1997). The objective of this 
paper is to summarize the research results from wind within 
a single tree to wind over a forest stand based on the 
studies of wind and management of a coastal pine forest in 
Niigata, Japan. 

Wind speed within a single tree crown 

Prediction of wind speed in a single tree crown as a 
function of position within the crown, crown structure and 
external wind is important for solving the problems of tree 
or stand damage by wind. The movement of air at all levels 
of a crown is important in the study of transfer processes 
among leaves, branches and their immediate environment 
(Murai et al. 1992; Zhu et al. 2000). Zhu et al. (2000) es¬ 
tablished a model of airflow within a single tree crown, with 
studies of Japanese black pine (Pinus thunbergii Pari.). 
Two Japanese black pine trees were examined and wind 
speeds at three vertical sections: VS<ns), oblique to wind 
direction; VS(pu, perpendicular to wind direction and; VS(pa>, 
along or parallel to wind direction) of the trees were meas¬ 
ured. 

In a conifer tree crown, the exponential form given in 
equation (2) is appropriate to describe the vertical wind 
profile (Zhu etal. 2000). 

U^^U^expi-a^-j-)] ( 2 ) 

/7 

where U m ( Z ), U ou t(z) are wind speeds within and without a 
tree crown at height z (m • s' 1 ), respectively, z is the interest 
height in the crown (m), H is the height of the tree crown 
(m), a s is a constant (dimensionless). 

Regressions were made to estimate the parameter a s by 
the wind speed data observed from the vertical sections of 
trees (Zhu et al. 2000). For estimating a s accurately, it is 
desirable to measure more individual vertical sections. 
However, it should be noted that accurate measurement of 
Uin(z) is always likely to be difficult in practice, because an 
instrument mounted in a tree crown is subject to the effect 
of local variability in the tree crown. The presented results 


indicated no significant difference among the vertical wind 
profiles, which were obtained from the total average, center 
average and the 3 vertical sections, respectively (p <0.05). 
Therefore, the wind measurement only at the center range 
(0.5 m around the stem) can provide the value of a s accu¬ 
rately for an individual tree. The measurement should be 
arranged in perpendicular to, or along with the wind direc¬ 
tion (Zhu etal. 2000). 

Therefore, wind speed within a single tree crown at 
height z (Ui n < 2 )) was estimated from the vertical wind profile 
by giving a constant value of a Sl together with a measure¬ 
ment of Uout(z). However, the wind speed Ujn(z) estimated 
from the vertical wind profile should only represent the 
mean wind speed at height z of the trunk, leastwise the 
center theoretically. It is necessary in some cases to esti¬ 
mate the wind speed at other positions within the crown but 
not at the trunk. For this reason, a horizontal wind profile 
expressed using crown thickness (CT) is developed (equa¬ 
tion 3). CT as an indicator of tree crown characteristics, was 
defined as a horizontal distance from the points for wind 
speed measurement in the observed section to the outer 
edge of a crown in windward direction, which was meas¬ 
ured corresponding to the points of wind speed observa¬ 
tions at each height. 

U inU) /U oulizr = e 0 cxp(eC T ) (3) 

where £ is a parameter to determine wind speed within a 
single tree crown according to crown thickness (dimen¬ 
sionless); Eo is a constant, given based on the mean values 
of relative wind speed from the observation of the three 
vertical sections (dimensionless), it is different according to 
various heights (Lower, Middle and Upper) (Fig. 1). 

Trees with longer crown thickness indicated a more wind 
reduction within the crown. The differential ratio dU(in Z /dCT 
can be described from an empirical regression expressed 
as in equation (3). Therefore, the wind speed at a horizontal 
plane within the crown can be estimated by providing pa¬ 
rameter e and CT. 

The vertical and horizontal wind profiles within a single 
tree crown of Japanese black pine can be orderly described 
in terms of height z and crown thickness as exponential 
functions, respectively. If the assumption that the estimated 
wind speed within the crown from vertical wind profile rep¬ 
resenting the mean wind speed at height z of the trunk is 
reasonable, then parameter £ in equation (3) can be de¬ 
rived from a single crown thickness measurement at the 
height corresponding to the estimated wind speed (Uin(z))- 
This is because the wind direction can be observed from 
instrument, and the crown thickness corresponding to the 
height z (distance from the trunk to the outer edge of the 
crown along the wind direction) can be measured. Using 
the equation (3) with known parameter e, wind speed at any 
position within the crown at a certain height (Lower, Middle 
or Upper) can be estimated reasonably. 
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Crown thickness (CT) /m 


Fig. 1 Relationship between crown thickness (CT) and rela¬ 
tive wind speed for total average from the three observed 
sections (Three parts, Lower, Middle and Upper were divided). 


Wind profiles inside the forest stand 

The analytical model for air movement within plant 
canopies developed in the literatures (Plate 1971; Bergen 
1971; Landsberg and James 1971; Kondo and Akashi 1976; 
Cionco 1985) equates the canopy to a two-dimensional 
distribution of infinitesimal momentum sinks corresponding 
to the canopy surfaces (Bergen 1971). The average flow in 
a horizontal plane is determined by a balance between the 
average divergence of shearing stress at that level and the 
drag per unit volume exerted by the canopy surfaces. The 
modeling of airflow within plant canopies described by 
Landsberg and James (1971), Thom (1971), Kondo and 


Akashi (1976), Cionco (1985) and Amiro (1990a) as 

^ = ZV / 2 (z<H) (4) 

where U m (z) is wind speed within canopies (m-s' 1 ), H is 
height of canopy top and the nominal tree height (m), and 
D 0 is the drag-vegetative attenuation factor defined as 

D 0 =±pC,A (z) (5) 

where Cd is drag coefficient (non-dimension), A< Z ) is the 
derivative of the leaf index with height (non-dimension). 

The model of airflow within canopies was formulated as 
exponential function of height in equation (6). 

u mz, = u h exp[-cr(l - —)] (6) 

ti 

where Uh is wind speed (m • s' 1 ) at height of canopy top H; a 
is a constant coefficient related to momentum absorption 
(non-dimension). The canopy flow index (attenuation coef¬ 
ficient) a determines the form of wind profile within canopy, 
and defined as equation (7) by Landsberg and James 
(1971); as equation (8) by Kondo and Akashi (1976) and 
Cionco (1985). 


\ c d U„L 

V 6 K 


m 


(7) 


where L is the average foliage density related to the leaf 
index. 


a = 



( 8 ) 


where Uo is mixing length within the canopy (m). 

Studies show that for the wind profiles in the trunk layer, 
at least on some occasions, bulge (Oliver 1975) or secon¬ 
dary wind speed maximum (Gardiner 1994) is found in the 
trunk layer (Oliver 1975; Zhu eta/. 1998). This finding has 
sometimes been explained by some blow-through phe¬ 
nomenon (Landsberg and James 1971) associated with 
edge effects, large gaps of forest or non-horizontal site and 
so on. 

In this study, wind profiles within the canopies of thinned 
and unthinned coastal pine stands were plotted according 
to various wind directions (Fig. 2). 

Obviously, the influence of wind direction on wind profile 
was stronger in the trunk layer than within the canopy layer. 
In the trunk layer, greater changes of wind direction, rang¬ 
ing from 46 to 180 degrees true north, were found for both 
thinned and unthinned stands. There was evidence that a 
certain amount of secondary maximum wind speed existed 
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in almost all of wind directions for the thinned stand, namely, 
wind velocity actually increased below the height under 
canopy and caused a negative gradient (du/dz<0). How¬ 
ever, there was no obvious secondary maximum occurring 
in the trunk layer of the unthinned stand. The phenomenon 
may be attributed to the fact that there was vacant space in 
the canopy of the thinned stand through which pressure 
pulses or gust of airflow could permeate to the bottom and 
spread from there. 


a leeward direction. The mean value of the attenuation 
coefficient a was determined as 2.2 and 3.4 for thinned and 
unthinned stands, respectively. This means that the foliage 
density in unthinned stand is larger than that of thinned 
stand (Peltola 1996b). In the trunk layer, the wind profiles 
can not be expressed using equation (6) because of the 
negative gradient of wind speed. 

Wind profiles above the canopy 



Uin(z)/UH 


Range of wind direction 


0-45 

136-180 

271-315 


—K—46-90 

-181-225 

—B— 316-360 


—*—91-135 
—*— 226-270 
—0—Total 


The earlier numerical models used concepts of eddy dif- 
fusivity, mixing length (L 0 ) and the shearing stress ( t 0 ) to 
describe the wind profiles above the canopy (Taketa 1964; 
Thom 1971, 1972; Kondo and Akashi 1976; Cao 1983; 
Cionco 1985; Kaimal and Finnigan 1994; Asai 1996; Liu 
and Kotoda 1996), i.e., logarithm law of velocity. For an 
idealized uniform forest expanse, the wind regime is gen¬ 
erally assumed to have a uniform flow with a constant ver¬ 
tical shearing stress above the canopy in adiabatic condi¬ 
tions (Munn 1966; Bergen 1971). The shearing stress of 
momentum flux can be expressed as 


T u 



(9) 


where K m is the exchange coefficient of momentum 
(m 2 -s' 1 ), p is air density (kg-rn' 3 ), U is the wind velocity 
(m-s' 1 ), zis height (m). 



Fig. 2 Wind profiles plotted according to wind direction for 
coastal pine thinned (A) and unthinned (B) stands. 

In the crown layer, the difference between the thinned 
and unthinned stands in wind profiles influenced by wind 
direction could be compared quantitatively using the at¬ 
tenuation coefficient a (equation 6). The parameter a was 
smaller for the winds with the windward direction (ranging 
from 271 to 45 degree true north) than those for the winds 
with other directions, and it was unstable for the winds with 


K m =L u u, (10) 

where Lo is mixing length (m), u • is friction velocity (m-s' 1 ). 
They are defined as 

L 0 = kz (11) 


where k is von Karman’s constant (= 0.40), and 



( 12 ) 


Substituting the above expressions for K m in equation (10) 

(13) 


gives, 

dU u 


dz KZ 

The above mentioned assumptions imply that the 
shearing stress is constant above the canopy or a surface 
layer, and z is the only length scale determining the shear¬ 
ing stress. Integrating equation (13) gives the classical 
logarithm law of velocity of the constant-flux (equation 14, 
the same as equation 1). 


U 


u, 

K 


[Inf 


z~d 0 


)J 


(14) 
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The parameters u-, cfo and zq in the equation of logarithm 
velocity law are called as wind profile parameters. Of which, 
cb and z 0 can be calculated from equation (14) as 


d 0 = z-z 0 exp(-^) 

(15) 

KU t 


z-d„ 

(16) 

© 

II 

exp(——) 





Strictly, equation (14) is just fitted for neutral condition 
(Aloysius 1995). When it is used in other conditions, it must 
be modified to conditions with non-neutral stability as 
equation (1). 

Generally, for obtaining the reliable estimation of the pa¬ 
rameters of do and Zo from the logarithm velocity law, it is 
necessary to make the measurement at least above height 
of 30 zo. This is because the turbulent flow in the logarithm 
velocity law can be expressed as 


IL 

u. 


-In (j-) + B 
K k s 


(17) 


where k s is the equivalent roughness parameter (m), z is 
height (m), fiis a constant. 

The logarithm velocity law can be applied only above the 
height of k s . Coefficient Sin equation (17) was determined 
as 8.48 in a situation of turbulent airflow (Xie 1987). 
Therefore, the effective height of k s above the canopy of the 
forest can determined according to equations (14) and (17) 
as 

k s = exp(Bxz 0 ) = 30z„ (18) 

The determination of parameter u- was measured actu¬ 
ally with the sonic anemometers. Here u* is expressed as 



where u' is longitudinal velocity fluctuation (m-s' 1 ), and w'is 
vertical velocity fluctuation (m • s' 1 ). 

An order-of-magnitude of stability length (L 1 ), calculated 
from the sonic anemometer measurements, showed a 
neutral condition for the observation. Therefore, d 0 and z 0 in 
the equation of logarithm velocity law for the coastal pine 
forest were estimated from equations (15) and (16) using 
the friction velocity u-, respectively. A total of sixty 
10-minute average profiles were used to calculate the pa¬ 
rameters of do and z 0 . The estimated values for thinned and 
unthinned stands were listed in Table 1. There was consid¬ 
erable variation in the estimation of both displacement 
height (do) and roughness length (z 0 ). The displacement 
height d 0 varied between 7.7 m and 0.6 m for the thinned 
stand, and between 8.9 m and 1.3 m for the unthinned 
stand. The roughness length z 0 varied between 2.39 m and 
0.18 m for thinned stand, and between 2.06 m and 0.06 m 
for unthinned stand. This result illustrated that estimating 
the parameters of d 0 and z 0 from wind speed profiles just 


above the canopy of the coastal pine forest was difficult, 
especially with the limited observations. The most reliable 
estimates of do and z 0 in this experiment probably came 
from those observations during July 17 to July 31 of 1999 
when the airflow blew perpendicular to the coastal pine 
forest. The average values of do and z 0 were 5.74 m, 0.52 
m and 5.56 m, 0.23 m for thinned and unthinned stands, 
which gave the values of ab/H=0.72 and 0.79, zo/H=0.06 
and 0.03 for each of the stand, respectively. The mean 
values of do/H and zo/H are similar to those obtained by 
other researchers (Landberg and James 1971; Gardiner 
1994). 


Table 1. Wind profile parameters above the canopy of thinned 
and unthinned stands. 



Thinned stand 

Unthinned stand 

Period 

u* 

do ' 

Zo 

u* 

do 

Zo 

Observed durina Julv 

< 17 to Julv 31 of 1999 




Maximum 

0.57 

7.73 

1.08 

0.61 

8.88 

0.61 

Minimum 

0.25 

3.48 

0.18 

0.18 

1.30 

0.06 

Mean 

0.42 

5.74 

0.52 

0.35 

5.56 

0.23 

Observed durina November 1 to December 27 of 1999 


Maximum 

0.45 

7.00 

2.39 

0.44 

8.19 

2.06 

Minimum 

0.37 

0.64 

0.26 

0.32 

4.02 

0.15 

Mean 

0.39 

3.17 

1.26 

0.39 

6.63 

0.67 

Total Mean 

0.41 

4.46 

0.89 

0.37 

6.09 

0.45 


Wind speed at a forest area 

The determination of wind regime in a potential or exist¬ 
ing forest area is an important aspect of any site assess¬ 
ment. For predicting annual wind speed for a forest area, 
Hannah etal. (1995) made an attempt to relate wind speed 
to topographic and geographic characteristics of the site in 
order to avoid the necessity for long-term on-site wind 
measurements. They concluded that the major variables 
relating to wind speed were altitude, ‘topex’ (a measure of 
the exposure of a site), roughness length and the distance 
from the coast. 

Over a period of several decades, the extreme wind 
speeds can be analyzed using extreme value statistical 
methods such as Gumbel (type I extreme value), Frechet 
(type II extreme value) and reverse Weibull (type III ex¬ 
treme largest values) distributions (Gross etal. 1994; Simiu 
and Heckert 1996; Galambos and Macri 1999). The prob¬ 
ability of occurrence of any value of hourly-mean wind 
speed can be described by a Weibull distribution fitted to 
the entire record of hourly values (Milne 1992; Quine 2000; 
Zhu etal. 2001b). 

Wind data (sample interval was 5 s, 10-minute mean 
wind speed, wind direction, maximum wind speed and the 
corresponding direction) were collected above the coastal 
pine forest during the period from November 16 of 1998 to 
July 16 of 2000 without break (87 696 sets of wind data). 
The frequency of occurrence of 10-minute mean wind 
speeds greater than 2.0 m-s' 1 was examined using 
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two-parameter Weibull function (equation 20) with the 
10-minute wind data (Fig. 3). 

^ exp|| (20) 

c c c 

where P w is the probability of occurrence of wind speed 
greater than mean wind speed U (m-s' 1 ), an exponent b 
(dimensionless) and a constant c ( m-s' 1 ) are parameters 
of Weibull distribution. The parameters b and ccan be used 
to derive the parameters (mode and dispersion) of extreme 
value distribution (Simiu and Filliben 1980, 1982; Simiu et 
al. 1985). 

The fitting data were limited to the wind speed greater 
than 5 m-s' 1 to minimize the effect of the poorer accuracy 
of the 10-minute means at low wind speeds. The distribu¬ 
tion was fitted to have b=1.36 and c=2.56 m-s' 1 . 


than value x in any one year, Ug is the mode of the distri¬ 
bution (i.e., the most likely value of x) and 1/A is the dis¬ 
persion, i.e., a measure of the width of the distribution of x 
(Milne 1992; Simiu and Heckert 1996; Quine 2000). 

The extreme wind speed, U ex tr corresponding to a recur¬ 
rence interval or return period, T ex tr (in years) is given as 


In equation (22), it has been found that equating x to the 
square of wind speed (U e xtr 2 ) rather than wind speed (Uextr) 
gives a better fit to the wind speed maximum (Milne 1992; 
Quine 2000). Thus, integrate equations (21)-(23), the ex¬ 
treme wind speed (U ex tr) corresponding to an observation 
period of T ex tr years can be calculated from equation (24). 



Wind speed /m-s' 1 

Fig. 3 Probability of occurrence for different 10-minute mean 
wind speeds over a coastal pine forest. The measured data and 
the Weibull distribution fitted to these measurements are com¬ 
pared. 

Generally, a data set conforming to the Weibull distribu¬ 
tion will have a cumulative distribution of its annual extreme 
values (Galambos and March 1999; Quine 2000; Zhu et al. 
2001b). The distribution of extreme values, which is usually 
used for building design purpose, is described by a 
Fisher-Tippett Type I distribution of the form in equation 
(21) (Miller 1982; Grigoriu 1984; Cheng and Chiu 1985; 
Gusella 1991; Milne 1992; Galambos and March 1999; 
Quine 2000). 

P c (x) = exp[- exp(-y)] (21) 

where y is a reduced variate given as 

y = ~{ X - u G ) ( 22 ) 

P G (x) is the probability that an extreme value will be less 


UL=U 0 -\ lnf-tad-J-)] (24) 

4 l ex,r 

where the parameters Ug and A can be estimated according 
to the method using Weibull parameters b and c, which is 
suggested by Milne (1992), that is, (Ug) 1/2 /c is a simple 
function of b, and AUg is approximately constant for an area. 
Therefore, the function of (Ug) 1/2 /c is derived by a third 
order polynomial of b (equation 25) as suggested by Quine 
(2000), and AU G « 5.0. 



= -0.5903Z> 3 + 4.4345/7 2 - 11.8633^ + 13.5690 


c 


(25) 


Generally, over a period of m days, data were collected 
for the highest wind speeds, while these observations are 
not independent. A method is suggested by Simiu and 
Heckert (1996) and examined by Galambos and March 
(1999) for modifying the data to achieve stochastic inde¬ 
pendence. According to Simiu and Heckert’s method, the 
relationships between the Gumbel distribution parameters 
(Ug and A) and the expected value E(u), and standard de¬ 
viation a can be expressed as 

Ju7 = E(u)~ 0.57722 (—)o ( 26 ) 

K 


A = (^-)a ( 27 ) 

K 

The parameters of the Gumbel distribution were esti¬ 
mated in this way. Applying this function to the recorded 
data over the coastal forest provided an estimate of Ug 1/2 = 
11.8(m-s' 1 )and 1/A= 27.7. From equation (24), the extreme 
10-minute mean wind speed at the measurement height for 
a 50-year return period was calculated to be 15.7 m-s 1 . 
The result showed that the approximation involved in both 
methods was in reasonable agreement. The probability that 
the extreme will be up to 25.0 m ■ s' 1 in any one year period 
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was calculated from equation (21), and presented in Fig. 4. 
This result suggested that it was a reasonable estimate for 
the coastal forest area to use the method for engineering. 



Fig. 4 Probability of extreme 10-minute mean wind speed 
over the coastal pine forest, being in the range from 5 m s' 1 to 
25 m s' 1 for one-year and fifty-year periods. 

Summary 

Generally, wind profiles can be divided into two or three 
parts for a forest stand, i.e., within canopy (crown layer and 
trunk layer) and above canopy. There are many studies 
related to wind profiles within and above the various types 
of forest vegetations. This paper summarized the gener¬ 
alization of wind profiles in a coastal pine forest area. First, 
wind profile within a single tree crown was modeled through 
observations of wind speed within the single tree crown of 
Japanese black pine, and the wind speed within any posi¬ 
tion of the single tree crown could be estimated by a con¬ 
stant, crown thickness and wind speed outside the single 
tree crown. Second, wind profiles inside the coastal pine 
forest stands were established by exponential function with 
the attenuation coefficient. Meanwhile, the difference of 
wind profiles inside thinned and unthinned stands was ex¬ 
pressed using the attenuation coefficient as well. Third, 
wind profile parameters (friction velocity, roughness length 
and displacement) over the coastal pine forest canopy were 
determined according to the ultrasonic observation. At last, 
the extreme wind speed at the coastal pine forest area was 
estimated with the methods developed for building design 
purposes. By summarizing the wind profiles in and over the 
forest stands, we found that some aspects such as rela¬ 
tionships between wind profiles and the tree or stand 
characteristics, the simple methods for predicting wind 
profile parameters, and the applications of wind profile in 
studies of tree physiology, forest ecology and management, 
and the detail ecological effects of wind on tree growth 
should be studied further. Because wind has a very wide 
variety of ecological effects, and has played a very impor¬ 


tant role in the development of forestry. 
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